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In supersymmetric models, a correlation exists between the structure of the Higgs sector quartic 
potential and the coupling of the lightest CP-even Higgs to fermions and gauge bosons. We exploit 
this connection to relate the observed value of the Higgs mass ruh ~ 125 GeV to the magnitude 
of its couplings. We analyze different scenarios ranging from the MSSM with heavy stops to more 
natural models with additional non-decoupling D-term/F-term contributions. A comparison with 
the most recent LHC data, allows to extract bounds on the heavy Higgs boson masses, competitive 
with bounds from direct searches. 



I. MOTIVATION 

The quest for SUSY has taken an unexpected turn with the Higgs discovery at 125 GeV [1]. Indeed, it 
is well known that the supersymmetric contribution to the Higgs mass is at most {rn*f[^'^)'^ < m^, implying 
that a large portion of the Higgs mass Arnf^ > 86^ GeV^ must originate from symmetry breaking effects. 
Within the MSSM, for large stop masses, top/stop loops provide this necessary contribution, but only at the 
expense of naturalness, as the large loop effects needed to increase the Higgs mass also destabilize the EW 
scale. Experiments are therefore telling us that, if SUSY exists, it is either tuned, or it doesn't fulfill Occam's 
principle and that more complicated models, with additional contributions to the Higgs quartic, have to be 
considered. 

Still, a common feature of most SUSY models^, is the Higgs sector, containing at least a particular version 
of a two Higgs doublets model (2HDM). Mixings in this extended Higgs sector, lead to modified tree-level 
couplings between the lightest CP-even Higgs and the SM gauge bosons and fermions, and provides a dis- 
tinctive signature of SUSY, complementary to direct searches. While the latter remain the most favorable 
strategy for SUSY searches (in particular in the most natural SUSY realizations, where states associated with 
the stabilization of the electroweak (EW) scale are expected to be light), modified couplings could be the 
strongest evidence for SUSY in particular regions of parameter space, such as those with compressed spectra. 

Interestingly, in 2HDMs, a correlation exists between the Higgs mass and its tree-level couplings to SM 
fields. Indeed, any contribution to the Higgs quartic potential, necessary in SUSY models to increase the 
Higgs mass from its tree-level value up to the observed value of approximately 125 GeV, also changes the 
relation between mass and hypercharge eigenstates and modifies the couplings of the lightest CP-even Higgs. 
In this article we investigate this correlation in detail, showing how different models that accommodate the 
observed Higgs mass also modify their Higgs couplings. We then confront these expectations with the most 
recent LHC data [8]-[12], which we use to extract limits on the parameter space of such theories (in particular 
on rriA and tan/3). 

We first show, with a simple and intuitive analytical approximation, how Higgs mass and couplings are 
correlated in SUSY models (or 2HDMs in general) (section II). Then we study, in turn, the MSSM with heavy 
stops (section HI), the MSSM with extra non-decoupling D-terms (section IV), and the F-term contributions of 
NMSSM-like models (section V) , where we also discuss a general class of models beyond the MSSM (BMSSM) . 
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^ An exception is the model of Ref. [2] where the Higgs is the neutrino superpartncr and there are no extra Higgs doublets. 
Ref. [3] also proposes a model with one doublet only while in Refs. [1, 5] additional doublets have been studied. 
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In section VI we comment on how these conclusions are modified in the presence of sizable loop-effects due 
to light SUSY partners and we leave for Appendix I the details related to our global fits and for Appendix II 
a summary of the formulas used in our plots. 

II. THE HIGGS MASS/COUPLINGS CONNECTION 

Supersymmetry requires the existence of two Higgs doublets, iJi,2 giving mass to leptons and down-type/up- 
type quarks. Limiting our discussion to the third family fermions, which have the strongest couplings to the 
Higgs sector, we consider 

C D -YbHiqb - YtH2qt - YrHjr . (1) 

Only a linear combination of Hi and H2 obtains a vacuum expectation value (vev) v = 174 GeV; its couplings 
to SM fermions and vectors equal those of a SM Higgs. Any quartic contribution to the scalar potential for 
Hi and H2 introduces, in general, a misalignment between this linear combination and the mass eigenstates: 
this misalignment is responsible for a modification in the Higgs couplings. The best way to see this is in the 
basis h, H, where only one state (ft) has a vcv. The angle (3 denotes the angle between these states and the 
neutral CP-even components of the gauge eigenstates Hi,H2: 

h\ = cos I3h + sin /3H (2) 
/12 = sin Ph — cos PH. 

In this basis, the couplings Eq. (1) of h and H to fermions are, 

- cos PYb{h + tan PH)bb, - sin jSYtih - cot PH)tt, (3) 

where couplings to charged leptons have the same form as for down-type quarks. Now, consider a general 
contribution to the quartic of the Higgs potential written in terms of h, H, 

AV{Hi,H2) = +Sxh^ + 5h^H + S2h^H^ + hhH^^ + S^H^, (4) 

where the J's are given dimensionless couplings. The first term contributes to the lightest CP-even Higgs 
mass as 

Am^ = mxv^; (5) 

in order to account for the observed value m\ w 125 GeV, 

Am^ = mf «2 - m|(cos 2/3)^ > (86 GeV)^ (6) 

is needed. Interestingly, the same physics that is responsible for also generates a mixing between h and 
if, via the term 6h''H . that leads to a modification of the Higgs couplings, as illustrated in FIG. 1. We can 
quantify these modifications, in the limit where H is heavy, by integrating out the heaviest eigenstate from 
the relevant part of the Lagrangian 

C^-SK'H- J2 Yff-fH-"^H\ (7) 

f=t,b,T 

where Yj^ , the coupling of H to fermion f = t,b, r, can be read from Eq. (3). For large m,H we can solve the 
equations of motion of H, giving H ft! Sh^Yp ff /m'jj, and obtain the effective interaction 



(8) 
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FIG. 1: The mixing between h and H , induced by the quartic interaction Sh^H , modifies the couplings of h to the 
fermions w.r.t to its SM value. 

Now, notice that the equations of motion for H imply a small vev {H) w 2^/25 [v^ /rn^j), so that the expression 
for the fermion mass is modified accordingly and we can write the coupling of the physical Higgs h = h — ^/2v, 
normalized with its SM value y^^^ — mf/v, as 

cf = ^^^—, V-^~ 1-4(5^^. (9) 

^ mf/v Yi^-2Y"5^ Y^ ml, ^' 

' ' t J mfj J " 

Using Eq. (3) to read Y^'^ , we finally obtain 

Cfc,r « 1 - 4tan/3(5X- , 
ct «l + 4cot/3J^. (10) 

This simple, yet important, expression summarizes the goal of this work: any new physics that is responsible 
for the large Higgs mass Eq. (5) also affects the Higgs couplings to fermions. This approximate formula allows 
us to understand qualitatively how this connection works and predicts whether a given contribution to the 
Higgs mass results in an increase or decrease of the couplings to tops and bottoms/taus (similar methods have 
been used in Refs. [13-15] to study Higgs couplings modifications). Nevertheless, notice that in our plots we 
always use the exact expressions listed in Appendix H, rather than Eq. (10). 

Deviations in the Higgs couplings to vectors can be studied in a similar way, giving 

which is generally suppressed w.r.t. deviations in the couplings to fermions (we have checked that in the region 
preferred by data this statement holds at better then the 2 % level and deviations in cy can be ignored). 

In principle, complete analyses of Higgs couplings in a SUSY context should take into account possible 
modifications of the tree-level couplings to up-type quarks, to down-type quarks (and leptons) and to vectors; 
at the loop level extra contributions from light SUSY partners to the couplings to giuons and photons could be 
present, and in total generality also the possibility of an invisible decay width should be considered (see Ref. [2] 
for a motivated scenario were the Higgs can decay invisibly in a SUSY context): a total of six parameters (see 
Refs. [16, 17] for a list of recent analyses of this type). Nevertheless, ignoring the last possibility, Eq. (11) 
tells us that in the simplest SUSY models, couplings to vectors are not expected to deviate much from the 
SM ones (this is not true when the Higgs sector is extended to include extra states in different SU{2)l 
representations that can mix with the Higgs, as we shall discuss in section V A) . Furthermore, the null results 
of direct SUSY searches suggest that SUSY partners should have masses of a few hundreds GeV and that 
their loop contributions to the effective hgg and /177 couplings might be small (we comment about this in 
section VI). For these reasons, in what follows, we orient our analysis mostly to the Higgs couplings to tops 
and to bottoms/taus and compare theoretical expectations with data through an intuitive simplified scenario 
where only ct,Cf, are free to vary, and all other couplings are fixed to their SM values. 
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III. THE MINIMAL SUPERSYMMETRIC STANDARD MODEL 

The technique of the previous section can be apphed also to the tree-level contribution of the Minimal 
Supersymmetric Standard Model (MSSM)^. The only contribution to the quartic potential comes from the 
D-term which, for the SU{2)l x U{1)y MSSM gauge group, reads 

AVmssm = = '-^^ {{c} - slfh' + 8(4 - sl)s,cph^H + ■■■) (12) 

with cp = cos/3 and sp = sin/3 and in what follows we shall also use tp = tan/3. This defines 

= ^(4-4f^ (13) 

2 

S - "^spcpicl-sj). (14) 

From Eq. (10), this gives 

2 

1 - — 4-sin4;9tan/3 (15) 

2 

1 + sin 4/3 cot /3. (16) 

which coincides with the usual decoupling limit of the MSSM [18] with the identification mn ~ ttia (which 
is accurate for mA,H S> mz or in the large tan/3 limit), and we will use in what follows in the comparison 
between exact and approximate results. At the same time, Eq. (13) provides the well known contribution 
to the Higgs mass = m\ cos^ 2/3; this tree-level result is modified by loop effects, in particular from top 
quarks/squarks, which we consider in what follows. 




A. Top Squarks with no mixing 

We begin with the case of top squarks with no mixing (realized in popular SUSY breaking mechanisms such 
as gauge mediation and gaugino mediation where a small trilinear coupling is expected [19]). The dominant 
loop contribution to the scalar effective potential is [18, 20], 

/:^Vstop^^\H2\\ (17) 



where. 



^log[mimijM^] (18) 



(a more accurate expression can be found in Appendix II). After rotating into the basis of Eq. (2) one identifies 

-^A = 4y (19) 

6 = -As^cp^. (20) 

From Eq. (18) and from Eq. (19) it follows that, in order to obtain a Higgs mass compatible with experiment, 
multi-TeV stop masses are required. Such heavy stops also destabilize the EW scale through loop effects and 



^ In this case, h and H can be thought of as the eigenstate of the mass matrix before electroweak symmetry breaking. 
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FIG. 2: Theoretical expectation for Higgs couplings deviations for the MSSM with heavy stops and no mixing, taking 
rrih = 125 GeV, showing contours of constant rriA (solid blue) and tan/3 (dashed), obtained from the exact expressions 
of Eqs. (68,69) of Appendix II. Also shown are the 68% (green), 95%(yellow) and 99%(grey) C.L. regions obtained by 
a global fit of the most recent LHC Higgs data, as explained in Appendix I, neglecting loop contributions to the hgg and 
/i77 couplings. The dashed red lines show the approximate results of Eq. (21) for mn = 300, 500 GeV. 



push the MSSM into fine-tuning territory [21]. Ignoring for a moment this tension, we can assume these loop 
contributions to be uniquely responsible for the large value of the Higgs mass, and write the deviations of Cb^t 
induced by loop effects Eq. (20) together with the ones from the tree-level potential Eq. (14), as 

TO? — TO?, cos 2/3 

9 '71? — TO 7 cos 2/3 , ^, 

Ct ~ l-(cot/3)2^^ f ^. (21) 

This shows that, in the MSSM with no stops mixing and for tan/3 > 1, the deviations in Cb (ct) are always 
positive (negaitive), as already observed in Ref. [15]. For large tan/3 the deviations in ct are suppressed, while 

/154GeV\^ 

(cb - 1 « . 22 

V J 

We can compare these results with the exact ones of Fig. 2, which shows the intuitive {cb, Ct )-plane mentioned 
above, and compares these theoretical expectations with the most recent data [8]-[12], using the methods 
described in Appendix I. We assume a heavy sparticle spectrum, that does not affect the Higgs couplings 
to gluons and photons, other than through Eq. (21) (this is motivated by the fact that in this example, we 
are assuming multi-TeV stops). Masses mn ^ 250 GeV can be excluded, almost independently of tan/3, as 
suggested already by Eq. (22) for a sensitivity to the hbb coupling of about 50%. In Fig. 3 we also show 
the CMS bounds on the traditional MSSM m^,tan/3 plane (for a recent analysis see Ref. [22]) from direct 
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FIG. 3: Exclusion plot in the mA,tan/3 plane for the MSSM with heavy stops (red), for models with additional non- 
decouplmg D-terms (blue) and F-terms (green); regions to the left of the lines are excluded. The shaded region corre- 
sponds to bounds from direct searches [59]. Left: present data; right: longterm projection based on [39] assuming no 
deviations from the SM, shaded region from Ref. [40] (the dashed part of the line corresponds to a region where As is 
bigger than 2 and non reaches the non-perturbative regime below approximately 10 TeV [21, 23]). 



searches of the heavy Higgs decaying into r pairs, as performed by CMS [59]. As can be appreciated, analyses 
of the Ught Higgs couphngs offer a complementary search strategy in the intermediate tan/3 region. 



B. Top Squarks with mixing 

In the presence of sizable A-terms, L and R top squarks can mix, inducing additional contributions to the 
Higgs effective potential [20, 25], 

Ay™" = ^\H2\^ + (y |i/li?2|' + X7\H2\^HiH2 + c.c), (23) 

where the values of A2, A5 and A7 depend in particular on the parameter fx and the trilinear At and their 
expression, at the one loop level, can be found in Appendix H. In the point of 'maximal mixing', when the 
trilinear term is \At — /iCot/3| = ^/&7ni (where is the geometric mean of the lightest stop masses), the 
contribution to the Higgs mass proportional to A2 is maximized, while A7 = 0. Recasting the potential in the 
h, H basis gives. 




where it can be seen that for large tan/3 (which is necessary in the MSSM to maximize the tree-level mass), 
the dominant contribution to the Higgs mass still comes from the first term A2, similarly to the case with no 
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FIG. 4: Same as FIG. 2, but for near maximal mixing and, again, we adjust ^mi^ mi^ £ [550,2000] GeV in order to 
obtain the observed Higgs mass. We take Xt = \/6±0.1 for the blue/red curve m order to show the influence, for large 
tan/3, of small deviations from maximal mixing; fj, = 400 GeV. 



mixing discussed in the previous paragraph. As mentioned above, this term is maximized by large mixing, 
with drastic effects and the stop mass can be as low as 550 GeV in this case. Nevertheless, a fine-tuning at 
the percent level persists due to the fact that large At terms also contribute to the Higgs mass-parameter [21]. 

Unfortunately, for a generic choice of /i and At , the multitude of parameters introduced by mixing weakens 
the Higgs mass/coupling connection as shown by Eq. (25) where sizable A5_7 can affect the Higgs couplings 
without contributing to the Higgs mass. We show this effect in Fig. 4 where we consider small deviations from 
maximal mixing: departures from A7 = X^axMzx _ g ^^.^ enhanced at large tan /3 > 20 and the contribution 
to 6 and to our predictions can be seizable. Nevertheless such large values of tan /3 are already in tension 
with rare B processes, such as Bg — >■ lJ,'^fJ,~ [26], and with direct searches for H/A — >■ fr [59], so that we do 
not expect our results to change significantly in the intermediate tan/3 region, where our bounds are more 
competitive, see Fig. 3. 



IV. EXTRA D-TERMS 



As discussed above, a 125 GeV Higgs in the MSSM is generally associated with fine-tuning. This suggests 
that the principle of SUSY, if realized at low energy in a natural way, extends beyond the MSSM, with 
new tree-level effects contributing to the Higgs quartic. The first possibility is to envisage additional gauge 
symmetries that contribute to the Higgs quartic, similarly to the MSSM gauge group [19, 23, 27]. In this 
section we study the example of an additional abelian gauge group under which Hi and H2 have opposite 
charges (as compatible with the /i-term). Then, the extra contribution to the Higgs sector quartic"^ 

AV^n{\H'i\'-\H',\')' (26) 



^ The form of the potential in Eq. (26) holds also for the non-abelian extension considered in Refs [23, 27]. 
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where, 

-= 'K. ■ (27) 

Here is the soft SUSY breaking mass of the MSSM singlets that breaks the \](l)x group (with gauge 
couphng gx) and Mz' the SUSY-preserving mass of the gauge boson. Eq. (27) shows that, in the hmit 
Mz' ^ rn^, the Z' can be supersymmetrically integrated out and the D-term contribution of the U(l)x 
group decouples: non-decoupling D-terms require a large soft mass ^ Mz' and result in an effective hard 
breaking in the Higgs sector. 

The contributions to 5\ and 5 are similar to Eqs. (13,14), with the substitution m^/v^ — > 4k. In the absence 
of other effects that affect the Higgs mass (we assume the loop effects of Eqs. (20,24) to be subdominant), we 
can fix K in order to obtain the observed Higgs mass ^, we can then write 

c, .1 + 2^^ (28) 



1 

1' 



ct « 1-2^-^. (29) 

meaning that, for tan/3 > 1, positive (negative) deviations are expected in Cb (ct). For large tan/3 the 
modifications in ct vanish, as usual, while those on Cf, asymptote to Cf, — 1 w (176 GeV/rri/f )^- This is shown, 
using the exact expressions from Appendix H, in Fig. 5. Differently from Fig. 2, the global fit of Fig. 5 includes 
the effect of a light stop at 500 GeV (as opposed to the previous section, where heavy stops were necessary 
to increase the Higgs mass, here this is taken care by the additional D-terms, and the stops can be naturally 
light, see also Section VI). Masses niH ^ 300 GeV can already be excluded, with better results in the small 
tan/3 region (see also Fig. 3). 

In principle we could relax the assumption that Hi and H2 carry equal and opposite U{\)x charges. In this 
case, however, additional structure is needed in order to generate a /i-term. For example an extra SM singlet, 
charged under U{l)x can generate this term by aquiring a non- vanishing vev. This extension, however, implies 
additional contributions to the quartic potential from F-terms which, as we comment in the next-section, are 
expected to dominate. 



V. F-TERMS, THE NMSSM AND THE BMSSM 

It is tempting to parametrize these new effects using an effective field theory approach with an expansion 
in powers of the scale of physics beyond the MSSM (in the example of the previous section, this would be the 
mass of the new gauge bosons Mz')- The most general such parametrization, however, lacks any predictive 
power (peculiar directions in parameter space can be found where an increase in the Higgs quartic coupling 
doesn't imply modifications of the couplings [28]). Nevertheless, as shown in Ref. [29], the leading order effects 
in such an expansion have a very specific form^ : 

C.^j d^e{^^^{H,H,f+Z^^{H,H,f^ (30) 

where Z = O'^mgusY is a dimensionless spurion that parametrizes SUSY breaking. This leads to additional 
contributions to the scalar potential, 

AFs = 2eiHiH2{HlHi + ff^ffa) + (2{HiH2f + c.c (31) 



* Notice that as tan/3 — ^ 1, all contributions to the Higgs mass from D-terms vanish; hence these expressions have to be trusted 
only away from this singular point: in FIG. 5 we show curves of constant gx (in the limit of large 2> M^') to show that 
in the region of interest the parameters are under control. 

^ For large tan /3 interactions at higher order in the expansion could be enhanced and dominate. 
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FIG. 5: Higgs couplings deviations in the MSSM with additional non- decoupling D-terms to raise the Higgs mass to 
125 GeV (on top of the effect of light stops mi = 500 GeV). The global fit (Colors as in Fig. 2) includes the effect of 
a 500 GeV stop (to be compared with Fig. 2 where the effects of stops on the fit are varnshingly small). 

with ei = Ai/i*/M and €2 = —X2msusY /M . We obtain 

6x = fisin2/3+^sin2 2/3 
4 16 

6 = -—cos 2/3- —sin 4/3. (32) 
2 8 

By construction /3 G [0,7r/2] and for the first term to contribute positively to the Higgs mass, a positive 
ei is necessary, implying an enhancement of Cb and a decrease in ct, similarly to the case studied in the 
previous section. The term proportional to £21 on the other hand, reduces for tan/3 > 1 the hbb coupling 
while increasing the coupling to top quarks, oppositely to the effects of D-terms. This is an interesting case 
that corresponds to the non-decoupling F-term contribution of an extra singlet, interacting with the Higgs 
sector via the superpotential term W — \sSH1H2, as in the NMSSM. Indeed, in the limit where the mass of 
the singlet is large, its contribution is given by the second term of Eq. (30), where M = Ms {msusY = ms) 
is the supersymmetric (SUSY breaking) mass of the singlet, and A2 — A| (notice that the singlet also gives a 
generally subdominant contribution to the first term of Eq. (30) with ei = — /j,*A|/(2Ms), which we ignore 
for the time being). 

If the largeness of the Higgs mass is due to a combination of the MSSM D-terms effects of Eq. (12) and 
the present contribution from F-terms due to the singlet (i.e. with negligible contributions from loop-effects). 
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FIG. 6: Coupling deviations in the NMSSM assuming a Higgs mass of 125 GeV in the limit where the singlet is heavy 
and it doesn't mix with the Higgs, but its contributions do not decouple. Global fit as in Fig. 5. 



then the Higgs couphngs to fermions are modified as 

^ ''"^"r" (33) 

e*-l + 42 ''^^ (34) 

which, for large tan/3, gives deviations in the hit coupUng of order Act « (60 GeV/mn)^, and in the couphngs 
to bottom quarks Acf, « t^(60GeV/m/f)^. We show the exact couphng deviations in Fig. 6 (we assume, again, 
the presence of 500 GeV stops, see section VI) where we also emphasize curves of constant As: values below 
^ 0.7 are perturbative up to the GUT scale, while for values 0.7 ^ Ag < 2 the non-perturbative regime 
is reached above a scale of 10 TeV [21, 30]. The bounds on mn that can be extracted from this analysis are 
very much dependent on tan /3, as can be seen in Fig. 3. 

While the approach of Eq. (30) encompasses large classes of models, its applicability is limited to cases with 
widely separated scales, such as the NMSSM where the singlet has both a large SUSY preserving and SUSY 
breaking mass^. In the opposite case, however, its interactions with the Higgs sector can induce mixings with 
the lightest GP-even Higgs and the analysis changes completely, as we now discuss. 



Triplets with hypercharge y it 1 and superpotential W = \TTH2H2 + XfTHiHi have also been considered in the literature: 
in the non-decoupling limit, their contribution to the potential is 

AV = \\t\'^hI + \\j.\''Hf (35) 

and 

^^ = ^4 + ^4' S=\^t\'4s^ + M^s$c^. (36) 
For large tan/3 only the H2 term is important and the results coincide with those of section III A. 
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FIG. 7: Feynman diagrams illustrating modifications in the couplings between the Higgs and the fermions (or vectors) 
due to mixings with the singlet S. 



A. Doublet-singlet mixing 

When the singlet is not much heavier than the EW scale, the above analysis ceases to be valid; moreover 
singlet-Higgs mixing can affect our discussion of section II (see also Ref. [31] for other LHC bounds on this 
possibility). Indeed, in this case, the potential Eq. (4) includes in particular the term AV{Hi,H2,S) D 
Sgsh? /2. Once h gets a vev, this term leads to a mixing between h and 5* so that h becomes a linear 
combination of the three gauge eigenstatcs: Eq. (2) now must include 

h = cos 6l(sin ph^ - cos + sin Os. (37) 

The mixing 9 can be estimated by using the techniques of section II: the term {Ss/2)h'^s corrects the h- 
propagator when s is integrated out and two of the h legs are replaced by vevs, as illustrated in the first 
diagram of FIG. 7. This correction, beside modifying the quartic structure S and 6\ as discussed so far, it also 
universally affects all h couplings by shifting the kinetic term to 

C,ff^{l + 2^-^)l^,h^'^h. (38) 

III Q Zi 

Indeed, making this kinetic term canonical leads to a universal suppression of all h couplings by the factor 

r2 2 

cos0«l-^, (39) 
where 9 is defined by Eq. (37) and the couphng of Eqs. (10,11) become, 

Cb ~ l-4tan/3(5^ - (40) 
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ct « 1 + 4cot/3(5^ - (41) 

cv^l-'-^. (42) 

Notice that in this case, if the singlet is light or if its couplings to the Higgs sector are large, sizable modifi- 
cations of the hZZ and hWW vertices can be produced. In principle, it is still possible to exploit the Higgs 
mass/coupling connection to fix 5 and then a simultaneous measurement of cy and Cb^t would allow to extract 
information about mn and about the mixing with the singlet. In practice, however, models of this type intro- 
duce many new contributions to the Higgs quartic potential and the Higgs mass/coupling connection looses 
most of its predictive power. We show this in the example of the NMSSM [32], where the superpotential 
W = XSH1H2 + kS^ /3 generates the following relevant terms in the potential (which add to the usual MSSM 
D-terms Eq. (12)), 

AVnmssm{Hi,H2,S) D mlS^ + X' {\H^H2\'' + S^\H^\'' + S^\H2\'')-{\AxSH^H2+\nHiH2S*^ +h.c.) (43) 
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where we assume real coefficients for simplicity {mg,Ax are soft SUSY breaking terms [32]). After the singlet 
obtains a vev (S) = vs, we can integrate out its real part, with mass ms, and obtain the effective quartic 
potential 



(44) 

(where we have neglected higher order terms in the couplings) and the mixing term 

(M2 - Ail sin 2/3) 

5s = , (45) 

with ^1 = A^A + 2\k,vs, and 1x2 = 2A^z)s- In this procedure, also contributions from the second diagram of 
Fig. 7 are taken into account. As usual the quartic potential can be written in terms of h, H and we find, 

5, = ^sin^2/3-i4 (46) 
16 8 mi, 

X ■ Aa /^i<5sCOs2/3 

= — — sm 4p 2 ^ (4 ' ) 



m 



2V2 



As it could have been foreseen, the multitude of parameters that characterize this model breaks the con- 
nection between S and S\ and it becomes possible to raise the Higgs mass independently of a modification 
of its couplings. Even for small itih a conspiracy between the MSSM D-term and these additional F-terms 
could allow for a large Higgs mass without any observable effect in the Higgs couplings (a similar example in 
the context of D-terms is discussed in Ref. [28]). 

Nevertheless, perturbativity up to the GUT scale (up to 10 TeV) limits the size of A < 0.7(2) and the 
necessity of a positive contribution to the Higgs mass from Eq. (46) , imposes an upper bound on the negative 
contribution proportional to Sg/m'g, as we show in the left panel of Fig. 8. Since the latter governs the 
coupling modification due to mixing through Eqs. (40-42), we see that in the perturbative NMSSM only small 
deviations are expected due to mixing, Acy ^ 5% for ms ^ v. Deviations in the couplings c^^t are still 
proportional to the parameter S which, as mentioned above, is now independent of the Higgs mass and would 
allow only to constrain the ratio fii/mn, which is not particularly interesting. 

In ASUSY[33], on the other hand, deviations can easily be of order unity. In particular, if 5 in Eq. (47) is 
positive (notice that for tan/3 > 1, both sin 4/3 and cos 2/3 are negative) we have Ct > cy > Cb, which enhances 
the rate of both /i — > 77 and h — > VV. Notice that if we consider only deviations in the tree- level couplings, 
an enhancement of /i — >■ 77 only, would require cy ^ Ct ^ ct, a region which is not touched by this model. 



VI. LIGHT SUSY PARTNERS 



So far we have studied modifications of the direct couplings between Higgs and fermions, restricting our 
attention to the 2HDM structure of the Higgs sector. When comparing with data, however, the presence of 
light sparticles can introduce additional nuisances, as they contribute via loop-effects to the hgg and hjj 
effective vertices. Naturalness suggests that only the partners of third family fermions be light (other bounds 
on natural SUSY have been studied in Refs. [34-36]); these have also the strongest couplings to the Higgs 
sector and have potentially a bigger impact than other sparticles. While staus and sbottoms have a negligible 
effect, light stops can change the analysis considerably [14]. We show this in Fig. 9, where we compare 
99%C.L. contours, assuming that a stop has been found, with mass m^ = 160, SOOGeV (dotted, dashed), 
with the contours without taking this effect into account (solid) ^. Since Ct itself affects Higgs physics mostly 



Recall that in our plots of the MSSM Figs. 2-4, since heavy stops are needed to increase the Higgs mass, we have assumed 
heavy stops in the fit too, with no seizable loop contributions to hgg and ft.77; in Figs. 5-8, on the other hand, where a natural 
spectrum is allowed thanks to the D-term/F-term contributions to the Higgs mass, we assumed = 500 GeV. 
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through a modification of the Higgs-gluon effective vertex* , the leading effect of light stops, which themselves 
affect the hgg effective coupling, results in a shift along the direction of Ct- As it can be seen, stops heavier 
than about 500 GeV have negligible influence on the fit. Nevertheless, in Fig. 9 we also show the global fit 
treating the stop contribution as nuisance and marginalizing over it: this is useful to take into account the 
possibility that a very light stop lies in a region inaccessible to direct searches. 

The mass of charginos is also directly related to the EW scale if the chargino is mostly Higgsino: then , 
In principle light charginos introduce an additional unknown through their contribution to the /177 coupling. 
However, for this to have any impact, small tan/3 < 5 [14], very light charginos m^± ^ 250 GeV [37] and 
large wino-chargino mixing (which is typically suppressed by inverse powers of the wino mass m^/M| [38] 
) are necessary. We consider this a peculiar, rather than representative, point in parameter space and we 
assume these effects to be small in our analysis. 



VII. CONCLUSIONS AND OUTLOOK 



In the MSSM, the tree-level Higgs mass is too small to account for the observed value of « 125 GeV. 
We have shown, using a simple analytical method based on an expansion in inverse powers of the heavy 
Higgs mass, how the physics that contributes to increase the light Higgs mass, also modifies the couplings 
of the lightest CP-even Higgs with fermions and gauge bosons. In the simplest examples (MSSM with no 
mixing between top squarks, MSSM with extra non-decoupling D-terms/F-terms) this connection provides 
distinctive predictions for the Higgs couplings, which allow us to extract bounds on the parameters niA and 
tan/3, competitive with bounds from direct searches [59]. Deviations in the couplings hZZ and hWW are 



Ct enters also directly through a contribution of a few percent of the tth associated production channel to the total production 
crossection and through the exclusive pp — > htt — > bbtt channel, which is however badly measured at present. 
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Ct 

FIG. 9: 99% C.L. contours from a global fit to the parameters Ch and ct, taking into account the loop effects due to a 
stop quark with = 160, 500, GeV (dot- dashed, dashed); in the solid line this effect is not taken into account, while m 
the colored contours stops (color coding as in Fig. 2) are treated as a nmsance and their contribution is marginalized 
assuming > 160 GeV. 

expected to be small in 2HDMs and we could use the intuitive Cb,Ct plane to show our results. In this way 
we could extract bounds on the heavy Higgs mass mn, depending on the model and on the size of tan/3, as 
we summarized in Fig. 3. 

Models that include extra gauge singlets that mix with the Higgs sector, can in principle be studied in a 
similar way. In the most popular realizations, such as the NMSSM, however, the large number of parameters 
of the model weakens the Higgs mass/coupling connection and predictability is compromised. Yet, theoretical 
consistency of the model, based for instance on the requirement of perturbativity, can strongly constrain the 
size of the expected effects. 

As long as the uncertainty in measuring Higgs couplings is dominated by statistical errors, more data will 
lead to better measurments. In the long term, with an integrated luminosity ~ 300fb~^, the sensitivity to the 
Higgs couplings to bottom/top quarks is expected to reach 15% [39, 40] and, as we show in the right panel 
of Fig. 3, some deviations from the SM are expected if mu^A ^ 400 GeV. At the same time, direct searches 
would have probed a much larger region of parameter space, but the bounds from Higgs couplings will remain 
competitive in the intermediate tan /3 region (better results can be achieved by considering ratios of couplings 
[41]). 

Let us conclude with a comment regarding bounds from flavour physics. The cross-section for fi 
processes is proportional to tan^ (3 and therefore this measurement practically excludes the region tan (3 >?>Q 
(depending on other parameters of the model [26]) while it has a relatively small impact for intermediate and 
small tan/3; in this regime the bounds discussed in this paper can be considered complementary. Constraints 
from 6 — S7 can be more important [46], but a fair comparison is difficult, as the amplitude for this process 
depends as much on the details of the sparticle sector as it depends on the parameters of the Higgs sector, which 
we consider here (in regions where the former are small, bounds on Type II 2HDMs exclude mfj± < 300 GeV, 
independently of tan/3 [48]). In any case, while the b ^ bounds are competitive with the Higgs coupling 
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bounds at present®, the latter are expected to become stronger as the integrated luminosity increases. 



Note Added: While this work was being finalized (see [43]), Ref. [44] appeared, which also considers 
modifications of Higgs couplings in SUSY models with additional, non-decouplings F-terms or D-terms. 
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APPENDIX I: Details of experimental fit 

To perform the statistical analysis of the data we take the following reescaled couplings of the higgs to the 
SM particles 



Ct 



yt 



vi 



SM 



Cb 



Vb 



Cv 



9hVV 

9hvv 



(49) 



but fixing Cv — ^ the reason of which is explained before in the text. 

We also take into account the stop loop effects which appear in the gg h, h 77, h gg production 
and decay modes. The main contributions to these effects [55-57] are given by the following expressions: 



CT(gg ~> h) 



T{h ^ gg) 



f^igg ^ h)sM r(/i gg) 



SM 



j^gg ^ ^93 ^ j^gg 
*i '2 



Agg 



(50) 



+ 



(51) 



^{h ^ 77) 
T{h -> 77)sAf 



All I Afl 

^w,SM ^ ^t,SM 



(52) 



1.28 a -0.28 q 



4mr 



4mx 



4mr mi 



(53) 



® An interesting contribution to b — >■ 57, which we have neglected throughout this work, comes from loop-effects involving 
squark-charginos or sbottom-gluino, which do not decouple when the mass of the superpartners is large and induce a coupling 
of H2 to down-type quarks, yb{H^bb + ei,H2bb) that can be enhanced at large tan/3 (which is interesting for the MSSM and 
for the D-term case of section IV). This can strengthen (weaken) the bounds from b s-y for negative (positive) ej, [45]. At 
the same time, however, it can affect the hbb coupling [14], 



Act 



tan (9- 



(48) 



which, for the MSSM or for the D-terms of Eq. (28), goes also towards strengthening (weakening) our bounds. Thus in the 
MSSM (also with additional D-terms) stronger constraints from Higgs coupling data are correlated with stronger b —> s-y 
constraints and vice- versa [47]. 
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The statistical analysis is performed using the latest signal strenght data given by Tevatron and ATLAS, 
and the one given by CMS at ICHEP. We didn't take the latest CMS data due to the fact that only a 
combination of 7 and 8 TeV is given in the signal strenghts at a higgs mass different than the one of ICHEP 
which doesn't allow us to extract them separatedly 

The signal strengths are assumed to follow a Gaussian distribution and we fit the data by minimizing a 
with the theoretical prediction for the signal strenght: 

i _ Hp Ct, Cj,, Cr)Cp BRi{a, Ct,Cb,Cr) , 

In the few cases where the given signal strenght errors are not symmetric we symmetrize them in quadrature. 
Statistical and theoretical errors are summed in quadrature without taking into account possible correlations, 
this approach is reasonable since at the moment the effect of this correlations is still small and can be neglected. 
In the other hand when comparing our fits using ICHEP data for both CMS and ATLAS we find an agreement 
of better than %10 between our figure and the one provided by them (we find this agreement to become better 
depending on the assumed cuts in the channels where they are not completely specified). 

The data used can be found in table I where all the channels taken into account are specified. In this 
table we see that for each channel a particular set of cuts is defined. These refer to the values of found in 
equation 54 which are the cuts for each higgs production mode: gluon fusion (G), vector boson fusion (VBF), 
associated production with a vector boson (A) and associated it (tth). Expanded this can be seen as: 



/L^p ^ p Sp 



'^gCg + ''V BfCv Bp + rrvv hCw H'^'^ ZhCz H'^'^'t'tHCltH 



where the cut efficiencies Cp for each production mode p corresponding to each channel i are reported below 
in table II. 



APPENDIX II: Details of the exact theory computation 

The most general two Higgs doublet model (2HDM) potential for the neutral components of the doublet is, 

+ (y|i/li/2|' + A6|i?lpi/ii/2 + A7|il2|'i?li?2 + C.c) (55) 

We have used the convention of Ref. [20, 25]. We can rewrite this potential in the h-H basis to obtain, 
(5 - ^c%Sfi + ^slcp + ^(Ag + A4 + A5)s^c^(4 - c2) + ^(c|(4 - cD + 2slcl) + y (4(4 - cD ~ 2slcl) 
= y4 + y4 + J(A3 + A4 + X,)sl4 + ^{4s,) + ^{slc,) (56) 
We will now give the values of A1-A5 in the different models we have considered. In the MSSM we have, 

mr 



^'-2^ ^'-2^ ^'-'2^ (^^^ 



An alternative approach can be found in [17 
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TABLE I: CMS, ATLAS and Tevatron data for the most sensitive channels. The cuts are classified as inclusive (I), 
associated production (A), vector boson fusion (VBF) or else ("fyx), see Appendix for details. denote the best 

fits for the 1.96 TeV Tevatron, and the 7,8 TeV LHC data. 
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At the one loop level we get the following additional contributions to the effective potential from top squark 
loops, 

AA2 = J^(6(2 + C2i/.);. + (l + C2i/s)a;tatyt^(12-xtat)) 
AA,5 = -j|^(l + cnZ.)(A)'a;?y' 
AA7 = -/ij^^xt(6 - xtat)(l + cai^s) 



Here, 



(58) 



Is = log[Af^^/mt], ?/t = TOt(mf)/(usin/3), 

xt = {At - n cot 13) /Ms, fl=^fi/Ms, at At /Ms 

cii = ^il2yf-32gl{mt)) 

C21 = ^^(62/? - 325^(TOt)) 

C31 = -^(92/2- 3252 K)) (59) 
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TABLE II: Cut efficiencies for production modes [11, 12, 49] of the channels of table I. Numbers in brackets give the 
efficiencies at 8 TeV, the others at 7 TeV and the overall normalization in each line factorizes. 
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For the BMSSM we get, 

AAs = 2ei AAe = 2e2 AA7 = 2e2 (60) 

For the D-term extension we get, 

AAi ^ 2k AA2 = 2k AA3 = -2k (61) 

Finally for the NMSSM with no doublet singlet mixing we get, 

AA4 = |AsP (62) 

The case of NMSSM with doublet singlet mixing has been dealt with in great detail in Sec. V A. For F-terms 
from triplets we get, 

AAi^2|AtP AA2 = 2|Aj.p (63) 

We can now write the mass matrix elements of the CP-even sector in terms of these couplings, 

M12 ^ 2v^[{\3 + X4)si3Ci3 + XqcI + Xrs}] - m\si3Ci3 (64) 
M12 = 2v^[Xic} + 2XGCpsp + X5sl]+m\sl (65) 
M22 = 2v^[X2sl + 2X7Ci3Sf} + X5cl]+m\c} (66) 
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where we have used [25], 

"^12 = Si3Ci3{m\ + 2\vjv'^ + Xf^fp^v"^ + Xjiisv"^). (67) 

Once we know the CP-even matrix elements we can easily find the exact coupling deviations. First we demand 
that the light Higgs mass, 

1 



ml = - i^Mn + M22 - ^{M22 - Mn)^ - ^Ml^j (68) 

is equal to 125 GcV by choosing an appropriate value of the stop mass in the MSSM, an appropriate value 
of gx for the D-tcrm extension and an appropriate value of A for the NMSSM. Now we can compute Cb = 
— sin a/ cos j3 and q = cos a/ sin /?, where a is extracted from 

tan2a=— — — . 69) 

M11-M22 
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